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A single-time quantum transport equation, which includes effects beyond the 
Cf") ■ quasiparticle approximation, is derived for Fermi-systems in the framework 

of non-equilibrium real-time Green's functions theory. Ternary correlations 
f-~» \ are incorporated in the kinetic description via a cluster expansion for the 

self-energies (e.g., the transport vertex and the width) truncated at the level 
of three-body scattering amplitudes. A finite temperature/density formula- 
tion of the three-body problem is given. Corresponding three-body equations 
reduce to the well-known Faddeev equations in the vacuum limit. In equi- 
librium the equation of state contains virial corrections proportional to the 
third quantum virial coefficient. 



I. INTRODUCTION 



A reliable description of the dynamics of Fermi-systems is commonly achieved in the qua- 
siclassical (QC) limit, which is based on the concept of a classical particle moving in the 
mean field and along a classical trajectory between successive instantaneous collisions. The 
QC picture assumes a hierarchy of time-scales; the time-scales on which the renormalization 
of the excitation spectra and collisions occur are expected to be smaller than the mean 
time-scale between encounters of dressed quasiparticles. A similar assumption applies to 
the length-scales. As well known, such a separation is possible in the limiting cases of a 
classical dilute gas and a highly degenerate Fermi-liquid (for a review see, e.g., |J0). The 
dynamical information in quasiclassical approximation (QCA) is contained in the single 
particle distribution function; the higher order distributions (determined form the BGKKY 
hierarchy) appear as time-dependent functionals of single particle distribution (for a review 
see, e.g., ||). A second element of kinetic description of Fermi-systems is the quasiparticle 
approximation (QPA), which assumes a sharp functional dependence of the energies of the 
quasiparticles on their momenta. The QPA is sufficient to constrain the time-scale (but not 
the length-scale) hierarchy to the QC limit; it does not, however, appear to be a necessary 
condition. 
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The need of extension of kinetic description of Fermi-systems beyond QPA arises in a 
number of physical situations. An example is the excited nuclear matter produced in heavy- 
ion collision experiments at energies ~ 100 Mev/nucleon (for a review see, e.g., M). The 
expanding nuclear system traverses a broad region in the density-temperature plane and, 
in particular, the temperature could be of order the Fermi temperature. In such situations 
the spectral function of the single particle excitations acquires a finite width. If the half- 
width of the spectral function is small compared with the amplitude, the corrections to 
the QP picture are small, and the hierarchy of the time-scale still holds (in other words, a 
distinction can be made between the renormalization effect leading to the broading of the 
QP peak and the elementary collision events between the single particle excitations with 
renormalized spectra). The finite width of the spectral function introduces in the system 
correlations beyond the renormalization of the single particle spectrum, which lead to the 
so-called virial corrections, and in particular to formation of bound states (in nuclear colli- 
sions - deuterons 0-0], tritons, alpha particles, etc.) and off-mass- shell propagation (such 
as Landau-Pomeranchuk-Migdal effect in the particle radiation processes [§]). Situation, 
similar to the one described above, is encountered also in the astrophysical context; e.g. 
in the finite temperature isospin asymmetical nuclear matter in the supernova collapse and 
proto-neutron stars. 

The picture above suggests seeking a kinetic equation for Fermi-systems that accommo- 
dates the separation of the time and length-scales (and hence the QC approximation) while 
allowing for effects beyond the QPA. [Formally this means that the time-scales of varia- 
tion of the self-energy functions (providing the renormalization and transport vertex) are 
small compared to the evolution time-scales for the propagators related to the quasiparticle 
distribution function] . 

Initial studies of the virial corrections to the Boltzmann equation were carried out using 
the density matrix approach || ; progress in this direction is represented by Snyder's kinetic 
equation |]10|| , where Bogoliubov's ansatz of the weakening of the initial correlations is used to 
close the BGKKY hierarchy at the two-particle level. The partition of the virial corrections 
between the collision integral and the drift term is not unambiguous fllj], however the 
kinetic equation reproduces the correct form of the lowest order (second) virial coefficient 
|T0| . The formation of the bound states and three-body processes within the density matrix 
approach has been considered e.g. by McLennan |I2| , Klimontovich et al [|H| and Ropke 



and Schulz |7] , who derive coupled kinetic equations for a gas viewed as a reacting mixture 
of atoms and diatomic molecules or nucleons and their bound states (the authors do not 
consider the third virial coefficient, however). 

The transport theory, which is based on the real-time Green's functions formalism and 



is largely due to Kadanoff and Baym |L4| and Keldysh |15|] , provides an alternative to 
the density matrix approach, and allows systematic diagrammatic approximations to the 
propagators and self-energies. The reduction scheme for the double-time integro-differential 
Kadanoff-Baym equations for the Green's functions to a single-time Boltzmann-type kinetic 
uses the QCA, which, as discussed above, separates the time and length-scales associated 
with individual collisions from those characterizing the inter-collision dynamics, and the 
QPA, which assumes a delta-shaped spectral function. The form and further approxima- 
tions for the scattering-in and -out rates, represented by the Feynman diagrams of elemen- 
tary scattering processes, vary considerably depending on the specifics of the dynamics and 
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composition of the underlying system; cf., e.g. 

The quantum kinetics beyond QPA has been addressed in a number of recent 

works [|19| -|22||. The memory effects in the non-Markovian kinetic equations, as shown in 

20fl , provide a sufficiently general description of correlations in quantum systems; the sec- 



ond order virial corrections emerge when the time-retarded spectral function is expanded to 
the next-to-leading order in retardation. This form of second-order quantum virial correc- 
tion agrees with the equilibrium result derived for semiconductors ||23||, Coulomb plasmas 



23 1 and finite temperature nuclear matter P5 |. The second order virial corrections in the 
local-time description and for the electron-impurity system have been given in pT[ in terms 



of the next-to-leading order expansion of the spectral function with respect to the width 



22 



The solution of the three-body problem for an isolated system is well know 
represented by the Faddeev |26|] or AGS 



and is 
p7| equations. The Kadanoff-Baym formalism 



has been applied to derive kinetic equations including triple collisions by Bezzerides and 
DuBois p8| ; their study, however, has been restricted mainly to the Born approximation, 
and consequently the collision integrals contain singularities due to the non-Fredholm type 
kernels of the scattering amplitudes. A different motivation for studying the three-body 
equations comes from the problem of the imaginary nuclear potentials, where the 2-particle- 
hole channel dominates ground state correlations in nuclei p5[| . A self-consistent random 
phase approximation to the respective three-body problem at zero-temperature has been 
given in ref. |p9]] , while the the finite temperature discussion of the latter problem is 
contained in refs. |~7|,|32| as an example of the application of the cutting rules for the 
multi-contour ordered Green's functions due to Danielewicz WJ\. In the near-equilibrium 



situation the nucleon-deuteron system has been studied at finite temperatures in the AGS 
formulation in refs. 13BIJ33I . 



Despite the progress mentioned above, the subject of the kinetics with three-body corre- 
lations needs further attention since (i) one would like to arrive at a kinetic equation which 
in the equilibrium limit leads to third order quantum virial corrections to the equation of 
state, (ii) the three-body collision processes should be free of singularities, therefore the 
appropriate Faddeev or AGS amplitudes should be used in the collision integrals (see also 
|34j]); (iii) the scattering amplitudes in the medium should reduce to the three-body Faddeev 
equations which are exact in the vacuum limit. 

The purpose of this paper is to address these items and to provide a link between the 
kinetic description beyond the QPA and three-body correlations. We adopt the familiar 
methodology of next-to-leading order expansion of the spectral function with respect to the 
spectral width, however instead of the decomposition of Green's functions into pole and off- 
pole terms [p0|-|22| , we use a decomposition in powers of the width. In the introduction to 
Section 2 we discuss the double-time kinetic equation and the gradient approximation to this 
equation (readers familiar with the formalism may skip this part of the section; the notation 
is as in Kadanoff and Baym [Q). Further subsections introduce the approximations to the 
spectral function, the method of decomposition of the diagrams in the QP and correlated 
parts, the Kadanoff-Baym ansatz and, finally, the single time kinetic equations. Section 3 
derives the three-body equations at finite temperatures and densities. The transport vertex 
and the collision integrals are discussed in Section 4. The equilibrium limit of the kinetic 
regime and the third quantum virial coefficient are derived in Section 5. Section 6 contains 



a brief summary of our results. 

II. DYSON EQUATIONS 

Consider a non-relativistic Fermi-system interacting via two-body forces. The Hamilto- 
nian of the system reads 

H = ^ £ / d 4 xV^l{x)V^ a {x) 

+ -E/ d 4 xd A x , 'il}l(x)-4)l,(x , )V fftff i(x,x r )'il} ff >(x , )^ a (x), (1) 

where i\) a {x) are the Heisenberg field operators, x = (f,t) is the space-time four vector, a 
stands for the internal degrees of freedom (spin, isospin, etc.). We shall use the matrix form of 
arrangement of time-contour ordered real time Green's functions, which is particularly suited 
when initial correlations are absent. In this case the upper branch, c_, runs from — oo to oo 
and the lower branch, c+, from oo to — oo. The initial correlation are commonly incorporated 
via the imaginary-time piece of the contour which accounts for the Kubo-Martin-Schwinger 



boundary condition ||17|| . We shall neglect this piece in the further consideration. The 



single-particle propagator matrix is defined as 

= f-i(Tif) ffl (xi)ilfl 2 (x2)) i(ipt 2 (x 2 )i) ai (xi)) \ (2) 

- { -»(^(xiX(x 2 )> -t(fl> n (xi)tf n (x a )))' [) 

where the indexes 1, 2... collectively denote the variables (xi, <7i), (X2, 02)..., the averaging is 
over a nonequilibrium state of the system; T and T are chronological and anti-chronological 
time ordering operators, respectively. 

The perturbative expansion for the Green's functions G can be arranged, as well known, 
in a manner similar to the ground state Feynman diagrammatic expansion utilizing Wick's 
theorem. The single-particle propagator, therefore, obeys the Dyson equation 

G 12 = Gg } + G ( S ® S 43 ® G 32 = G { $ + G u ® S 43 ® Gfl (3) 

where the superscript (0) refers to the free propagator, ® stands for matrix multiplication 
along with the folding over the internal (repeated) indexes and the self-energy matrix S has a 
structure identical to eq. (^|). The integro-differential form of the Dyson equation is obtained 
by acting the evolution operator G,^ 1 = id/dt\ + V\/2m on eq. (^) (the energy scales are 
measured relative to the fermion chemical potential) and employing the free-particle Dyson 
equation, G^ 1 (g> G 12 = c 2 S^(xi — x 2 ); here and below ctj (i = x,y,z) are the Pauli 
matrices and the subscript refers to the evolution operator of free particle. Subtracting 
the resulting equation from its conjugate one finds 

( G Q2 l * - Gq! 1 ) ® G 12 = G 13 ® S 32 ® a z - a z ® S 13 ® G 32 . (4) 



The matrix structure of the Green's functions and the self-energies is not the optimal one yet. 
The retarded and advanced functions are preferable [fL5,17,28 , since they obey an integal 
equation in the lowest order of QCA. In terms of real parts of these functions (Re G R = 
Re G A = Re G and Re T, R = Re E" 4 = Re S) the kinetic equation, e.g., for the G < component 

is 



Gq3 — ReS 



13 



G 



32 



-[ReG 



13 



^32J 



1 



1^13' ^32/ 



1 



1^13' ^32/ ) 



(5) 



where [ , ] and { , } denote the commutator and anti-commutator, respectively, and a sum- 
mation (integration) over repeated indexes is assumed. If the dynamics of the system per- 
mits a quasiclassical treatment, i.e., if the characteristic inter-collision length-scales are 
much greater than the inverse momenta and the relaxation times are much larger than the 
inverse frequencies, QCA can be applied by separating the slowly varying center-of-mass 
four-coordinates from the rapidly varying relative coordinates. Performing a Fourier trans- 
form with respect to the relative coordinates and keeping the first order gradients in the 
slow variable the QC Kadanoff-Baym kinetic equation reads 



i{ReG~\ G<} pB +*{£<, ReG} PB = S < G > - £ 



>G<, 



(6) 



where all functions depend on the four energy-momentum vector p = (u, p), and the center 
of mass space-time vector x. Here the four- component Poisson bracket is defined as 



= df_ dg_ _ d£ dg_ _ d£ dg_ df_ dg_ 

Ouj dt dt Ouj dp Or dr dp 



(7) 



Equation (H) is the starting point of the reduction of double time kinetic equation to the 
single time form; one recognizes the left-hand-side as a generalized drift term, while the 
right-hand-side as the gain and loss terms of the collision integral. 

A. Spectral decomposition 



The Dyson equation for the retarded/advanced Green's functions to the leading order in 
QCA has an integral form, with the formal solution 



G R/A ( P ,. 



CO 



Z R/A (p,x)±ir] 



(8) 



where the free single particle spectrum is e p = p 2 /2m. The respective single particle spectral 



function, defined as a(p, x) = i G R (p 



,x 



a(p,x) 



G (p,x) , has the general form 



[uj-e p (x)} + [j(p,x)/2] 



2 ■ 



(9) 



where j(p, x) = i 
e p + ReS(p,x)| w=£ 



S K (p, x) — Ti A (p,x) 



is the width of the spectral function and e p (x) = 
is the dispersion (the pole value of the retarded Green's function (El)). 
Following refs. |23|-|25|, we expand the spectral function up to the next-to- leading order 
term in the power series expansion with respect to the width 7 



V 

a(p, x) ~ 2nz(p, x)S(uj — e p {x)) — *y(p, x) — — -, (10) 



UJ — Ep{X) 



with the short-hand notation V'/(uj — £ p (x)) = d/duj [V/(uj — e p (x))]. Note that the self- 
energy appearing in the denominator of the second term of eq. ( 110) via the dispersion relation 
is restricted, to first order in 7, to the mass-shell. The wave-function renormalization, z(p, x), 
in the same approximation is 

n , fdoj' , , , V 
z{p,x ) = 1+ /— 7 (w,P,x — , 11 

J 2tc uj' — uj ^=£p 

where we used the integro-differential form of the Kramers-Kronig relation: 

—-Rei:(u,p,x) = / ImZ(u ,p,x) = - / — j(u , p, x) -. 

OUJ J IT CO — UJ J Z7T UJ — UJ 

One may observe, combining expansion ([10|) and (|ll|), that the spectral sum rule, 

J—a(p,x) = l, (12) 

is fulfilled to any order in the expansion with respect to the width. 

B. Correlation bracket 

The reduction of the double-time kinetic description to a tractable single time descrip- 
tion requires projecting out different energy states, yet keeping the maximal information on 
the spectral properties of the system. The QPA projects out just a single energy value (de- 
fined by the on-mass-shell condition) leaving out the off- mass- shell dynamics of the system. 
On the other hand, keeping the full spectral function in the kinetic scheme is prohibitive 
for numerical applications, and is not necessary if the off-mass-shell effects appear as cor- 
rections. The process of projecting out the energy states can be systemized in terms of 
simple decomposition rules (see also [[23|-|2~5|| ) . Consider an auxiliary step of decomposition 
of tu-integrated arbitrary off-mass-shell function(s) with respect to the spectral width. A 
substitution of the spectral function in the integral allows us to establish simple rules; e.g. 
for a function of a single energy argument one finds 

J a(uj)F(uj)duj = F(e p ) + {F(e p )} cmv , 

{^ P )} corr = /^7(w,P,x) J^\F^) ~ F(e p )], (13) 

where "corr" stands for correlation. A product of two functions of a single energy argument 
decomposes as 

J a{u)F x (u) F 2 (uj)duj = F 1 (e p ) F 2 (e p ) + F 1 (e p ) {F 2 (e p )} coii 

+ F 2 (e p ) {Fi(£ p )} corr + {^(eA F 2 (e p )} cott , (14) 



/tit ) "P' 

^7(w,P,x) -^— [FjM - Fi(e P )][F 2 M - F 2 (£ P )]. (15) 

For a single function of two energy arguments one finds 
r 1 

J a(Wl)a(W2)F(tJi,W2)dWldW2 = ^(%l,£p 2 ) + 7j ({ F ( £ Pl^P2)} CO rr + { F ( £ P2^Pi)} CO rr) > ( 16 ) 

— -n/(u u p,x) [F(u x ,e p% ) - F{e px ,£ n )], (17) 

and so on. Any diagram now can be decomposed in terms of these rules in the leading and 
next-to-leading order terms in the spectral width. The coherent terms of the type {Fi-F^jcorr 
assume one and the same deviation of the functions i*\ and F 2 from the mass-shell, and are of 
a higher order in 7 in the incoherent limit {Ti} corr {-^jcorr- The decomposition rules above 
imply expansion in orders of 7 and differ from those used in pl] , p2| where the decomposition 
is in the pole and off-pole terms. 

C. Examples 

Let us illustrate the use of correlation brackets by writing the two-body T-matrix and 
one-loop polarization in the next to leading order approximation. These functions, as well 
known, contain the full information on the two-body scattering in the particle-particle and 
particle-hole channels respectively; apart from the fact that they represent one of the few 
examples where the resummation series can be summed-up in a closed form, their significance 
is due to the fact that the pole of the T-matrix signals the onset of the superfluid phase 
transition (Cooper phenomenon), e.g. |T^Ji7| , |37|j38| , while that of the polarization function 



the onset of the growth of density fluctuations and liquid-gas phase transition, e.g. J39| , [40 
The contour ordered T-matrix equation reads: 



1234 — V 1234 



+ ^Vl278® G 75 G 86 ®T 563 4, (18) 



where V1234 = cr z V1234, is the time-local bare interaction. In general, each element of the 
matrix T is a function of four time arguments. The time locality of the potential, however, 
implies a double-time structure identical to eq. @. For the same reason, the particle-particle 
propagator product G G = Gw, p ] should be considered as a single matrix. The components 
of the scattering amplitudes, needed for complete specification of the self-energies, can be 
chosen as the retarded/advanced ones; the remaining components are provided by the optical 
theorem. For the retarded/advanced T-matrix the Bethe-Salpeter equation in the Wigner 
(or mixed) representation reads 

(p\T R / A (P,x)\p') = (p\V\p')+ i Y,(p\V\ P ")Gl / p f(p",P,x)(p"\T R / A (P,x)\p'), (19) 

p" 

where the leading order terms in the gradient expansion of the product g R / A T R / A have 
been kept. Here the subscript [pp] indicates the particle-particle channel and p, P rela- 
tive momentum and total four-momentum respectively. The two-particle Green's function 
appearing in the kernel of equation (|T9|) is given by 



Ggf (pi, Pi, x) = Y,H { G> ( p i/ 2 + Pi. x) G> (P 2 /2 - Pl , x) 

0)1 P 2 

-G < (P 2 /2 + p 1 ,x)G < (P 2 /2-p 1 ,x)}R R / A (P 1 ,P 2 ,x), 



where 



R R/A (PuP2,- 



(2tt) 3 S 3 (P 1 -P 2 ) 



(20) 



(21) 



Ei(x) — E 2 (x) ± it] 
is the free two-particle resolvent. The full off-mass-shell kernel now can be decomposed as 

G 2pf (Pi> p i> a;) - [l - /(£+, x) - /(£-, x) 



- {/(e+, x)} corr - {/(e_, x)} corr J R R / A (E 1 ,e + + e.,x) 
+ [1 - /(£+,£) - /(e_,x)] [{i?^( J B 1 , £+ ,x)} corr 

+ { J R i?/A (£ 1 ,£_,x)} co „ 

- [{/(£+, z)7e R ^( J B 1 , £+ ,x)} corr + {/(£_, x)^/ A (E,e_,x)} c 
+ O (^ 



(22) 

where abbreviations e + = ep 1 / 2+pi and £_ = £.Fi/2-pi has been used; the relation between 
the propagators G < , G > and the distribution functions / (i.e. the ansatz) is given in the 
next subsection by equations fl26|) and (p7|). The correlation brackets for the TZ function are 
more complicated than those defined above; here the convention is adopted to denote only 
the energy argument in the sum e + + e- with respect to which the correlation is constructed, 
i.e. 



{n R/A (E,e + )} c 
{f(e + )n R (E,e + )} c 



1 V 



,R/A 



K (&, oo + e_) - ^o ' (E, e + + e 



R/A. 



^ O) £+ — ^ 

[< /A (£, w + e_) - < /A (£, e + + e. 



(23) 



(24) 



The first term (in brackets) in eq. (f2"T| ) corresponds to two on-mass-shell propagating 
quasiparticles with intermediate state phase space occupied by both the quasiparticles and 
off-shell excitation. If the correlation terms are neglected one recovers the quasiparticle limit 
of the Bethe-Salpeter equation. Since in the equilibrium limit the correlated states obey Bose 
statistics one may conclude that, apart from the Pauli-blocking due to the quasiparticles, 
the intermediate state propagation is Bose-enhanced by the off-shell excitations. In nonequi- 
librium situations, of course, the occupations of quasiparticle and correlated states are the 
solutions of the respective kinetic equations (see equations flS5|) and ( p6[ ) below). The sec- 
ond term in brackets corresponds to a two-particle propagation, where one of these is in an 
off-shell state; consistent with the next-to-leading order expansion, the intermediate state 
propagation is suppressed by quasiparticle Pauli-blocking. To the next-to-leading order in 
7 the third term vanishes in the incoherent limit; in the coherent limit only an equilibrium 
treatment is possible. The kinetic theory, in fact, does not provide an equation determin- 
ing the time evolution of the off-mass-shell distribution function. It is worthwhile to note 



8 



that the hole-hole propagation (the G > G > term) is included in the two-particle propagator. 
When this term is dropped the T-matrix equation in the zero temperature limit reduces to 
a Brueckner-Galitskii p5| , |36| type integral equation for a slightly non-perfect Fermi-gas. 

The decomposition for the one-loop polarization function requires a decomposition of 
the intermediate state particle-hole propagator 



G*{ A ( Pl , P u x) = J2 E { G> ( p 2/2 + Pl , x) G< (P 2 /2 - p t , x) 

m p 2 

- G< (P 2 /2 + Pl ,x) G> (P 2 /2 - p x ,x)} R R/A (P u P 2 ,x). (25) 

Since one can trace a complete analogy to the T-matrix case we shall skip further details. 



D. The ansatz 

The kinetic equation (|6|) is still incomplete; it should be supplemented by a relation 
between the functions G K and G > . A natural choice is the Kadanoff-Baym ansatz, 

-iG < (p,x) = a(p,x) f(p,x), (26) 

iG>(p,x) = a(p,x) [l-f(p,x)}. (27) 

The spectral function has already been defined as a(p, x) = i G R (p,x) — G A (p,x) = 
i [G > (p, x) — G < (p,x)}, therefore the ansatz replaces one of the propagators by the func- 
tion f(p, x) which has the meaning of a quasiparticle distribution function; it reduces to the 
Fermi distribution in the equilibrium limit. The non-diagonal elements of the self-energy 
matrix can be expressed, similarly, via the spectral width and a distribution function, 

iE < (p,x) =7(p,x)/'(p,x), (28) 

-iH > (p,x)= 1 (p,x) [l-f{p,x)}. (29) 

Though generally f'(p, x) ^ f(p, x), these functions, however, differ by terms of higher order 
in the gradient expansion than is need for the present discussion (see ||19| ). The reduction of 



the drift term in equation (|6|) can be facilitated by introducing auxiliary Green's functions, 

- iG<^ (p, x) = 2vr f(p, x) 8(u - e p (x)), G< n (p, x) = [1 - z(p, x)] G<^ (p, x) (30) 

G<(p,x) = -G< n (p,x) + Z<(p,x) ^—-, (31) 

uj- s p (x) 

and similarly for the G > (p,x) function with particle occupation replaced by the hole one, 
fip) ~ ¥ [1 — fip)]- The measure of the change of quasiparticle Green's function due to the 
wave-function renormalization is included in the Gf en propagator, where "ren" stands for 
renormalization. In terms of these functions the expansion of off-diagonal Green's functions 
with respect to the width 7 is 

G<(p,x) = G<^ (p,x) + G<(p,x). (32) 

If this decomposition is substituted in the kinetic equation @, the drift terms corresponding 
to the leading and next-to-leading order contributions decouple: 



{ReG-\ G 



< _ f< 

7^0 rcn 



{ReG-\G< + G< n } pB = -{£<,ReG} 



P.B. 
13. 



S > G<-S < G 



<r» 



P.B. 



(33) 
(34) 



The decoupling of eqs. (^) and ( |34"D derived in a slightly different manner in ref. [f22| . 
justifies the common practice of dropping the term {E'S Re G} p B from the Kadanoff-Baym 
equations when taking the QP limit. 



E. Single time equations 

Now we are in a position to derive the single time kinetic equations corresponding to 
equations ( p3| ) and fl34|). Integrating these equations over u one finds 



d de p d 
dt dp dr 



de p d 
dr dp 



f(p,r,t) = I(p,r,t), 



l + £^) {/(p ' r ' t)}corr = | /t(p ' r ' t) + |: /r(p ' r ' t) ' 



(35) 
(36) 



where the collision integrals are 

I(p,r,t) = J ^[E>(w,p,r,t) G<(u,p,r,t) -Z<(u,p,r,t) G>(u,p,r,t)], (37) 

(kC(iJ P ' 1 [X > (u,p,r,t)G<(u;',p,r,t)-Z<(u;',p,r,t)G>(u;,p,r,t)], (38) 



It(p,r,t) 
I r {p,r,t) 



(2tt) 2 u-lu 
duj duj' V 







- \G<(u;,p,r,t)— E < (^,p,r,t)-S > (^,p,r,t) 
(2n) 2 uj — ijj' \ dp 1 

+ S< (u, p, r , t) A [ G > ( w ', p , r , t) - Gf< (w', p, r, t)' 



(39) 



Equations (^) and ( p6|) couple the time evolutions of the distribution functions for the on- 
mass-shell propagating quasiparticles and off-shell excitations, which are described by the 
distribution functions /(p, r, t), and {/(p, r, t)} CO rr respectively The terms involving space- 
time derivatives in (^) on the right- and left-hand side of the equations balance each other 
separately; e.g. the correlated distribution function can be determied from I t by a direct 
integration (provided the initial conditions are known). The latter collision integral involves 
only the on-mass-shell singlet distribution functions, whose time-evolution is governed by the 
kinetic equation (35|). One may now compare the kinetic equations ( |35"D and fl36|) to those of 
refs. pO|-p^] . Ref. |20| gives a single kinetic equation for the Wigner distribution function. 
Summing eqs. ( p5f) and (|36|) we find a kinetic equation which is consistent with the eq. (49) 
of ref. 
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except the terms involving the space gradients of the correlated distribution 
function. This difference lies presumably in a different partition of the collision and mean- 
field terms (note that the both terms are governed by the same timescales). Present form 
seems to be more suited for establishing the conservation laws; e.g. the number density 
conservation is obtained automatically by taking the momentum integral, while ref. |21J 
drops the space gradients. Refs. [0,0] present two kinetic equations; compared to their 
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equation for the 
Wigner function, 
Wigner function, 
evaluated in the 
since it includes 
can judge, differ 
Summing eqs 
conservation 





singlet distribution function, which is defined as the pole values of the 

our eq. ( 35 ) differs by the fact that it includes only the 7 — > part of the 

The second kinetic equation ref. |22[ gives in the integral form, which is 

Born approximation. The content of our equation (|36"D is again different, 



Of 



the dynamics of all correlations oc 7. The collision integrals, as far as we 
considerably in their general form (cf. eq. (56) in ref. |22|j). 

and (PB|) and integrating over momentum we find the particle number 



(40) 



[n frcc (r, t) + n COTT (r, £)] + V r [jft- ee (r, t) + j co A r i *)] = °, 



where 



n, 



n free (r,£) 
r,t) 

jfvee(r,t) 
Jcorr \f ; tj 



d 3 p 

(2^p 

d 3 p 
J^rj 3 

d 3 p de p 



(2tt) 3 dp 

d 3 p de p 

(2tt) 3 ~dp 



f(p,r,t), 

{f(p,r,t)} corr , 
f(p,r,t), 
{f(p,r,t)}< 



(41) 
(42) 
(43) 
(44) 



and the terms O^ 2 ) have been dropped in the expression for the current of off-shell ex- 
citations, Jcoir(f,t); (symmetrization of the collision integrals in the usual manner shows 
that they vanish after integration over momentum). The separation of the Wigner distri- 
bution function f tot (p,r,t) in leading and next-to- leading order terms, as noted above, is 
not the only possible one. In fact, the separation in the pole and off-pole terms has been 



with 



£2|, 


where 


. 










/tot(p) = 


r du 

J 2ixi 


K 


>o(^ 


p)+G<(oj,pj_ 






fosip, 


r,t) 


-/ 




P)- 



Z(p)f(p) + fas(p) 



(w,p)]. 



(45) 



(46) 



Though formally equivalent, the latter partition does not fulfill the frequency sum rule 
at each order of the expansion and one has to sum up at the least first two terms. The 
partition in powers of 7 has the advantage that it does so at any order of the expansion and, 
in addition, it permits a simple interpretation of various terms because, close to equilibrium, 
the distribution function of correlated exciatations have a Bose-type spectrum. 

III. THE THREE-BODY PROBLEM 

We proceed to formulate the three-body problem at finite temperatures and densities. 
The resummation series for three-body scattering amplitudes have the general form, 
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r=v+v®g®v 
= v + v®g ®T = v + T®g ®v (47) 

where Go and Q are the free and full three-particle Green's functions, and V is the interaction 
(we use the operator form for notational simplicity; each operator, as in the two-particle 
case, is combined in a 2 x 2 matrix, with elements defined on the contour). In the case of 
pair-interactons, the net interaction is V = V12 + V23 + V13, where V a p is the interaction 
potential between particles a and (3; note that the potential matrices are time diagonal due 
to the time-locality of the potential. The kernel of eq. ([|7]) is not square integrable: the pair- 
potentials introduce delta-functions due to momentum conservation for the spectator non- 
interacting particle and the iteration series contain singular terms (e.g., of type VapQoVap to 
the lowest order in the interaction). The resulting ambiguity is eliminated by a rearrangment, 
originally due to Faddeev ||26|| , which plugs the singular term in the channel ^T^-matrices 
corresponding to the case when the interaction between the particles a,/3 and the third 
particle 7 is neglected. The total T'-matrix has a decomposition 

t = t (1) + t (2) + t (3) , (48) 



where 



T {a) = V M + VVy ® g <g> T. (49) 



and a(3 r y = 123, 231, 312. The channel (a(3) transition operators resum the successive 
iterations with driving term V a p, 

Tap = V a( 3 + Va/3 ® Go ® Tap] (50) 

they are directly related to the two-body effective interaction in the a/3-channel, e. g. two 
particle T-matrix, eq. flT8|). In terms of these functions the total three-body amplitude is 
completely determined via the three coupled integral equations 

r {a) = r M + 7> 7 ® g ® (r m + r (7) ) , (51) 



which represent non-singular Fredholm-type integral equations; their formal structure is 



identical to the Faddeev equations in the vacuum p6fl . To see the statistical nature of the 
system (finite temperature and density) which enters via the contour ordering of the opera- 
tors reflected in their matrix structure, let us write eq. (|5T|) in the momentum representation. 

E. g. for the retarded component of T^ 1 ) one finds 



<fc23,<?i|T*«(M^|*4,<7i> = (fc 2 3,9i|T 2 f(M',*0|fc 23 ,<?'i> 




&23 Hi , Ql 

2 4' 23 2 



T R W(t,t"',K)\k'? 3 ,qZ) 



/k 23 3q x q ± 

+ {— + —,-k 23 + Y 

x(k^q'^it'",t",K)\k'; 3 ,q^) 

x ^ _ M, _ fc '; 3 _ £\T*(f, t\ K)\k' 23 , q[) 

x k^ 3 dq' 2 n dt'"dk'[ 3 dqldt'\ (52) 
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where we used the time-locality of the interaction; (the obvious dependence of the func- 
tions on x has been dropped). Here the momentum space is spanned in terms of Jacobi 
coordinates, 



I, 1, , 2 

3 J 



K = p a +p + p 7 , k a = -{p a -p ), q y = -(p a + p ) - -p 7 . (53) 



The free three-particle Green's function has different types of factorizations depending on 
the particle-hole content of the three-body T-matrix. One may identify four types of direct 
factorizations 

Go lppp] (ti,t 2 )=9(t 1 -t 2 ) [G>(t 1 ,t 2 )G > (t 1 ,t 2 )G > {t 1 ,h) - G<(t 1 ,t 2 )G<(t 1 ,t 2 )G<(t 1 ,t 2 )], (54) 

Qo[ PP h]{tiM) = 0{t 1 -t 2 ) [G > (t ll t 2 )G > (t ll t 2 )G < (t 1 ,t 2 ) - G < (t ll t 2 )G > (t ll t 2 )G < (t u t 2 )} , (55) 

So[ P hh](ti,t 2 ) = e(t 1 -t 2 ) [G > (t ll t 2 )G < (t ll t 2 )G < (t 1 ,t 2 ) - G < (t ll t 2 )G > (t ll t 2 )G > (t u t 2 )} , (56) 

Go[hhh](ti,t2) = e{t 1 -t 2 ) [G < (t l ,t 2 )G<(t 1 ,t 2 )G < (t 1 ,t 2 )-G > (t 1 ,t 2 )G > (t 1 ,t 2 )G>(t 1 ,t 2 )}, (57) 

where subscripts [p..] and [h..] refer to particle and hole states respectively. The first line 
corresponds the 3-particle channel which is dominant in Fermi-systems interacting via short- 
range forces, and to which we shall restrict ourselves. The Faddeev decompositon, eq. (p2| ) . 
and five exchange diagrams for this case are shown in Fig. 1. 

The signs of the exchange diagrams are given by (—1)™, where n is the number of the 
line intersections in the diagram. The respective exchange terms for the remainder of the 
diagrams are obtained by interchanging the time-direction of lines. 

A Fourier transformation in eq. (j52|) with respect to the relative times leads to the 
3-particle T-matrix 

(fc^lT^JQIfcWi) = (k 23 ,q x \r R (K)\k' 2 „q' x ) 

+ /{(^-^,-fc 2 3-||T^ 2 )(K)|^^) 

T R ^{K)\k'^q'i)} 

Q(K',k'{,,q'J,)5(K-K ! ) 



2 4 2 



x 



Q - Q' + ir] 



/ i ^13 ^2 . // Q 2 



T£{K)\Va,4i) 



x du'{ 3 dk" 3 di4dq^dn'dK', (58) 

with 

Q(K, k 13 , q 2 ) = {G> (K/3 + k 13 + q 2 /2) G> (K/3 - k l3 + q 2 /2) G> (K/3 - q 2 ) 

- G< (K/3 + k 13 + q 2 /2) G< (K/3 - k 13 + q 2 /2) G< (K/3 - q 2 )}, (59) 

and four- vector notation K = (Q, K), k = (uj, k), q = (u, q). The momentum representation 
of equations for T^ and T^ follows from ( |58"D by permutations of (Greek) indexes 1 <-> 2 
and 1 <-> 3. 
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FIGURES 
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Fig. 1 

FIG. 1. a) The three-body amplitudes in the 3-particle channel after Faddeev rearrangement; only one 
of the terms of the type T( Q ) appearing on the right hand side is shown, b) The direct and five exchange three 
particle propagators. The diagrams for processes involving hole propagations are obtained by interchanging 
the direction of the lines. 

To the lowest order in the expansion of the off-diagonal single-particle Green's function 
we find the QP result 



{k^ qi \T R U{K)\k'^c( x ) = (fc M ,gi|T*(JO|fc , asi qf'i) 



&23 

2 



3gi 

4 ' 



*»-£ 



T R W(K)\k'l 3 ,q'£) 



+ ( , — &23 + 

2 4 2 



x 



n 



q ^T R W(K)\k f l 3 ,q»)}Q(K,k>l 3 ,qZ: 

~2^ i 2 "' _ ^13 ~ ~~3~ ^23 1-^)1923' 9l) ,. ,, , ,, 

dk 13 dq 2 , 



S K/3+k'{ 3 +%/2 € K/3-k'^+<%/2 € K/3-% + l7 l 



(60) 



with 



Q(K,k 13 ,q 2 



= l-f(K/3 + k 13 + q 2 /2)-f(K/3 
- / (K/3 + fe 13 + q 2 /2) f (K/3 - k 13 - 
-f(K/3 + k 13 + q 2 /2)f(K/3-q 2 ) 

-f(K/3-k 13 + q 2 /2)f(K/3-q 2 ). 



- fe 13 + q 2 /2) - f (K/3 - q 2 ) 
ft/2) 



(61) 



Equations (0) and (|BTf ) provide the Faddeev amplitudes appropriate at finite densities and 
temperatures. The medium effects enter these equations in two ways: (i) the renormalization 
of the single particle spectrum in the resolvent of eq. fl58|), and (ii) statistically occupied 
intermediate state propagation. The first line of eq. (|B"TD is the Pauli-blocking of intermediate 
propagation of three particles and is the dominant term in the dilute limit. The remainder of 
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the terms, in the same limit, tend to Bose-enhance the intermediate propagation, which can 
be seen from the equilibrium identity f {u)\) f (uji) = gij^i + ^2) [1 — /(^i) — f( u 2)], where 
/ and g are the Fermi and Bose distribution functions. Inclusion of the next to leading 
order terms in the powers of 7 can be constructed in a manner similar to the two-body case, 
however would lead to implicite inclusion of four body correlations, which are beyond of 
the scope of this paper. In the limit Q — > 1 these equations reduce to the original Faddeev 
equations (of course with the particle spectrum e = p 2 /2m). 

IV. SCATTERING INTEGRALS 

A. Transport Vertex 

The self-energies in the matrix formulation of nonequilibrium theory have a structure 
topologically identical to the equilibrium case. We shall adopt the cluster expansion for 
the self-energy matrix [0,0] and, consistent with the discussion above, truncate it at the 



three-body level: 

5^12 = T1234 <8> G43+ + / 7~165432 <8> G43+ G 65 +, (62) 

where + means an infinitesimal later time point. The non-diagonal elements of the matrix 
equation provide the collision rates on the right hand side of the kinetic equations. In the 
mixed representation one finds 

X < (p l ) = -i52((p l -p 2 )/2\T<(p 1 +p 2 ,)\(p l -p 2 )/2)G>(p 2 ,) 

Pi 

+ J2[( k 2s,Qi\^ {1)< (K)\k 23 , qi ) +(k ia ,q 2 \T®<(K)\k u ,q a ) 

P2,P3 

+ (k l2 ,q 3 \T^<(K)\k l2 ,q 3 }}G>(p 2 )G>(p 3 ). (63) 

The expression for E > follows by simultaneous interchange of < and > signs. The Jacobi 
momenta in the three-body T-matrix are related to the arguments of the self-energy and 
the Green's functions via the relations fl53|). Equation ( |63"D does not imply the usual form 
of the transition probability as a square of the transition amplitude times the statistical 
factors (i.e. Fermi Golden Rule). For the purpose of recovering this form one may use the 
non-equilibrium optical theorem 



x (2n) 4 5 4 (P-p 3 -p 4 ), (64) 



(Pi|T<(P)|p 2 ) = -i J2(Pi\T R (P)\^-^) G<(p 3 )G<( P4 )(^-^\T A (P)\p 2 ) 

P3P4 Z Z 



(k 23 , qi \r<(K)\k' 23 ,q[) = Y: E (k 2 3,q 1 \T^ R (K)\k 45 ,q 6 )G<(p 4 )G<(p 5 )G<(p 6 ) 

a,j3 P4,P5,P6 

x (k A5 ,q 6 \T^ A (K)\k' 23 ,^)(27r) 4 S 4 (K -p A -p 5 -p 6 ). (65) 



Combination of the equations ( |63|) ,( |6^) , and ( |65D yields the usual form of the quantum 
Boltzmann collision integrals with scattering probabilities defined as 
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(pi,p 2 |W|p 3 ,p 4 ) = ((pi-p2)/2\T R (P)\(p 3 -p 4 )/2)((p 3 -p i )/2\T A (P)\(p 1 -p 2 )/2) 
x <5 3 (p 1 +p 2 -p 3 -p 4 ) (66) 

(Pl,P2,P 3 |>V|p4,P5,P6)=E( fe 23,qi|r^ iJ (K)|fc 45 , q6 ) ^45,^1^)^(^)1^23,^) 

x (2tt) 3 5 3 ( Pl + p 2 + p 3 - Pi ~ P5 ~ Pe), (67) 

where it is understood that the Jacobi momenta in the three-body transition probability W 
should be expressed through the momenta Pi, i = 1...6, via the relations (|53|). 

The diagonal elements of the matrix equation (|62l) provide the solution of the integral 
Dyson equation and thus, the renormalization of the single particle spectrum. Apart from 
the details (which can be deduced along the lines of the previous discussion, see Section 1 
E), let us note that (i) in the next-to- leading order approximation to the spectral function 
the mean-field due to the excitations occupying off-mass-shell states contributes to the the 
real part of the self-energy; (ii) an additional contribution to the self-energy comes from the 
off-mass-shell T-matrices folded by the on-shell quasiparticle distribution - a contribution 
which is missing in the usual Brueckner theory; (iii) the three-body processes contribute to 
the energy shift in the single particle spectrum via the second term of the cluster expansion. 

B. Virial expansion of collision integrals 

The general form of the collision integral can be arranged in the form of virial cor- 
rections in much the same manner as the diagrams of scattering theory (Section 1 C). In 
principle, in the next-to-leading order approximation this collision integral can be separated 
into the quasiparticle and correlated contributions in two different ways, corresponding to 
the partition in the pole and off-pole terms and in the powers of 7. The latter partition 
is physically more appealing, since it involves the coupling between the kinetic equations 
for on-shell quasipaticles and off-shell excitations. This decomposition for the two-body 
scattering integral contains five terms 

I(p, x ) = j(°) (p, x ) + J (1) (p, x) + I {2) (p, x) + J (3) (p, ar) + J (4) (p, x) ; (68) 

a similar decomposition emerges for the three-body collision integral X. The first terms in 
the decomposition i.e. the scattering integrals for the quasiparticles, 

/(0) bi) = J2 (P^P2\W\p3,Pi)S(e pi +£ P2 -e P3 -£ PA ) 

P2P3P4 

x {/(fh) /(P 2 )[l - /(ps)][l - /(p 4 )] - [1 - /(Pi)][l " f(P2)]f(Ps) /(p 4 )} , (69) 
j(0) (Pi) = H (Pi,P2,P3\W\p 4 ,p 5 ,p 6 )5{e Pl +e P2 +e m -e Pi -e 5 -e 6 ) 

P2P3P4P5P6 

x {/(pi) /(ft) /(p 3 )[l - /(p 4 )][l - /(p 5 )][l - /(Pe)] 

- [1 - /(pO][l - /(p 2 )][l - /(p 3 )]/(p 4 ) /(Pfe) /(Pe)}, (TO) 

describe the scattering of two/three incoming uncorrelated quasiparticles into uncorrelated 
outgoing states. Note that the retarded amplitudes in the collision probabilities include, 
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apart from the usual Pauli-blocking due to the uncorrelated quasiparticles, the off-shell in- 
termediate propagation effects and the phase space occupation due to the off-shell excitations 
in the medium. 

The remaining collision integrals describe scattering processes where one of the incoming 
or outgoing excitations is in a correlated state. (The upper index labels the correlated parti- 
cle). Note, however, that the processes involving correlated states implicitly contain at least 
an extra quasiparticle. The terms of higher than zeroth order in the three particle collision 
integrals, therefore, involve processes with four or larger number of quasiparticles and can 
be ignored. As to the two-body case, the result for the second term of the decomposition 
reads 



-Pi, 



/(1) (Pi) = J2 (Pi^P2\W\p 3 ,Pi) [{f(e pi )} COTT S(s pi + e P2 - e P3 - e f 

P2P3P4 

+ f( £ Pi){$( £ Pi)}corr + {/(£pj %pi)}co 

X {/(ft) [1 - f( P3 )] [1 - /(p 4 )] + [1 - f(p 2 )] f(p 3 ) /(p 4 )} ; (71) 



where the functions {S(e pl )} COTT and {f(e Pl ) S(e Pl )} corr are defined in complete analogy to 
eqs. ([23|) and (p4[). The explicit expression for 1^ follows from 1^ by interchanging 
Pi ^ P2, while those for I^ 3 ' and B 4 ' from simultaneous interchanges p$ <-> p 2 , p 3 «-> Pi 
and Pi «-» pi, p 3 <-> P2, respectively. Consistent with keeping only the first order terms in 
7, the scattering amplitudes in the integrals /W, i = 1,2,3,4, are those including only the 
quasiparticle contributions. The first terms in the bracket in eq. (|7"T|) describes scattering of 
excitations, when one of these is in an initial or final correlated state. This term describes a 
real process with energy conservation in the scattering event and irreversible deformation of 
the wave functions of the scattering particles. The second term corresponds to the virtual 
scattering of quasiparticles with explicit nonconservation of energy and reversible deforma- 
tion of the wave functions; the third term is of the coherent nature discussed above (see Fig. 
2 a). 

Had we preferred the decomposition in the pole and off-pole terms, then the zeroth 
order term would have acquired an additional factor which is a product of the wave-function 
renormalizations of all four quasiparticles. For the next-to-leading order virial correction 
one finds 

I {1) (Pi)= J2 (Pi,P2\W\p 3 ,p 4: ) j(e P3 + e P4 -e P2 ) 

P2P3P4 £ Pl ' £ P2 £ P3 £ P4 

x {/(^P3+ £ P4 -e*a)/(P2) 

x [1 - /(Pa)] [1 - /(P4)] - [1 - f(e P3 + e P4 ~ ^ 2 )] 

x [l-/(p 2 )]/(p 3 )/(P4)}. (72) 

The fact that this collision integral corresponds to an effective three-body process can be 
visualized by inserting in the latter equation the explicit expression for the spectral width 
(in the quasiparticle approximation) 
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a) 




Fig. 2 

FIG. 2. a) The virial decomposition of the collision integral in the two body case; the open circles 
denote the T matrix in the quasiparticle approximation, while the filled ones include also the off-shell 
contributions, b) The modification of ncxt-to-leading order correction to the gain part of the collision 
integral after substituting the explicit expression for the width; all the final state particles are now on the 
mass-shell. 



7(Pi) = - E {pi,p2\W\p 3 ,p4)(2Tr)S(e pl +e P2 - e P3 - e P4 ) 

P2P3P4 

x {/(p 2 ) [1 - /(p 3 )] [1 - /(p 4 )] + [1 - /(pa)] /(p 3 ) /(p 4 )}, 



(73) 



which gives 

i (1 \pi)= E E (pi,P2\w\ P3 ,p,) 



v 



P2P3P4 P5 PC P7 



£pj "T £p 2 ^P3 ^P4 



'P7V 



x |((pi-p 5 )/2|T (pi+p 5 ,£ P3 + e P4 + ^5 -epa,») (P6-p?)/2 / 

x (2vr) 4 5(p 3 +P4+P5-P2-P6- Pi) S(e P[i + e P4 + e P5 - e P2 - e P6 - e, 

x {/(^ 3 + £ P4 " £ P2 ) /(P2) [1 - /(Pa)] [1 - /(p 4 )] 

- [1 - /(£„ + e P4 - e P2 )} [1 - f(p 2 )] /(p 3 ) /(p 4 )} 

x {/(p 5 ) [1 - /(p 6 )] [1 - f(Pr)] + [1 - /(Ps)] /(pb) f(pr)}. (74) 



This result is essentially self-explanatory; the underlying process is illustrated in Fig. 2 b. 

Let us turn to the collision term in the kinetic equation (|36|) . The first term in the virial 
expansion is beyond the quasiparticle approximation, and e.g. for the collision term I t has 
the form 



ir(p)= E — — z — ^(p,p2\w\ P3 , Pa ) 

P2P3P4 £ "T" £ P2 E P3 S V4. 

x {/(p) /(p 2 )[l - f(Ps)][l - /(p 4 )] + [1 - /(P)][l - /(Pa)]/(Ps) /&*)}■ (75) 



Compared with the first term of the virial expansion of the regular collision integral I(p, x) 
this term differs by the off-shell energy denominator; in contrast to the virial correction 
terms here all four particles are in the final on-shell states, however their energies are not 
matched by the energy conservation condition. The latter expression defines the correlated 
part of the distribution function as a functional of the on-shell QP distribution function via 
integral form of the kinetic equation 



V. VIRIAL CORRECTIONS IN EQUILIBRIUM 

When integrated over momentum the kinetic equation provides number density conser- 
vation; in the equilibrium limit the conserved density is the sum of quasiparticle density and 
the density of correlated excitations 

ntotO) = ^frcc(^) + ^corrO) = J 77^ [/(%, v) + {f(Sp, r)} corr } . (76) 

Consistent with our approximations, the equation of state contains virial corrections up to 
the third order. Before giving the expression for the third quantum virial coefficient let us 
briefly recapitulate the second virial coefficient (see also 



A. The Second Virial Coefficient 

The starting expression is that for the width at the two-body T-matrix level 
7 (p0 = -2i£{((Pi -p 2 )/2\lmT R ( Pl +p 2 )\( Pl - p 2 )/2) G<(p 2 ) 

P2 

- (( Pl -p 2 )/2\T<( Pl +p 2 )\( Pl - P2 )/2)lmG R (p 2 )}. (77) 

First we simplify this expression by eliminating the T < function using the optical theorem 
(0) in equilibrium. The latter can be transformed, using the operator relation 2 ImT^ = 
i(T < — T > ), to the following form: 

( Pl \lmT R (P)\p 2 ) = -~ £ o(ws) aK)( Pl |r fi (P)|(p 3 - Pa)/2) [1 - f(u 3 ) - /(w 4 )] 

P3P4 

x ((p 3 - Pi )/2\T A (P)\p 2 )(2n) 4 S(P-p 3 -p i ), (78) 
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where f(u>) = [exp (P(u — //)) + 1]~ is the Fermi distribution function, p being the chemical 
potential, and f3 the inverse temperature. Introducing the distribution functions in the eq. 
(|64|), noting that in equilibrium 

fMfM = gfai + w 2 ) [i - /(wi) - fM], 

where g{ui) = [exp (j3(oj — //)) — 1]~ is the Bose distribution function, and comparing eqs. 
(||) and (|78D one finds that 

( Pl \T<(P)\p 2 ) = -2ig{u l + u 2 ) ( Pl |ImT R (P)|p 2 ). (79) 

The expression for the width j(p) simplifies to 

7(p0 = 2$>(a; 2 ) [^ + co 2 ) + f{u 2 )} ({ Pl - p 2 )/2\ImT R { Pl + p 2 )\{ Pl - p 2 )/2). (80) 

Substituting the expression for j(p), eq. (|T3|), in the expression for the correlated density 
and using the equilibrium identity 

[g(E) + f(cu 2 )] [f(E - u 2 ) - /(<*)] = ^(E) - g(ui + w 2 )] [1 - fM - /(w 2 )] , (81) 

one finds the following form of the two-particle correlation function |23|-f25|j 

•dE 



{f H (e Pl )} co „ = f^-g(E)b(e Pl ,E), (82) 



where 



b(e Pl ,E) = 2£[1 - /( %1 ) - /( %2 )] {((p! -p 2 )/2|ImT*(£)|( Pl -p 2 )/2) ^ReR R (E) 

P2 «^ 

+ Im^(E) -^(( Pl -p 2 )/2\lmT R (E)\( Pl -p 2 )/2)} (83) 

is the two-particle spectral function and R is the two-particle resolvent. 



B. The Third Virial Coefficient 

To derive the third virial coefficient we first relate the T <,> matrices to the imaginary 
part of the retarded T R - matrix. The optical theorem reads 

(fc 23 g 1 |Imr^(^)|fe45a 6 ) = \Y. £ (k 2 3q 1 \T R ^(K)\k 45 q (i }[G<(p 4 )G<(p 5 )G<(p 6 ) 

a/3 P4,P5,P6 

-G > (p i )G > (p 5 )G > {p 6 )](k 45 q 6 \T A ^\K)\k 23 q 1 ) 
x S±(K- Pi -p 5 -p 6 ). (84) 

where the momenta pi i = 1 . . .6 are related to k a p, g 7 via relations fl53|). Using the equilib- 
rium identity 
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fMfMfM = /(^l + uj 2 + w 3 ) [l - /(o;i) - /(w 2 ) - /(w 3 ) + fMfM 

+ f(u 1 )fM + fMf(u 3 )], (85) 

one finds in the QP limit: 

(fc 23 qr 1 |X < (ir)|fc 45 g 6 ) = 2i/ ( Wl + cj 2 + cu 3 ) (k 23 q 1 \ImT R (K)\k 45 q 6 ), 
(knqill^Wlknqe) =2i[l-f (^ + u 2 + w 3 )] (fc239i|ImT R (ir)|fc 4 5g 6 >. (86) 

The three-particle correlation function takes the form 
{/ m (^i)}corr = -2i /-^ X! (fc23qfi|Imr fl ( J K')|fe 2 3g 1 ) [/ (wi + e P2 + e P3 ) + g(e P2 + e P3 )} 

x [1 - f(e P2 ) - f(e P3 )} lf(u) - f(e Pl )] -?— (87) 

£ p\ ~ u 

Next we use the identity, valid in equilibrium, 

[/(w + e P2 + e P3 ) + g(e P2 + e P3 )} [1 - f{e P2 ) - f{e P3 )} [/(w) - /(e^)] 

= [/(#) - fie* + £ P 2 + £ P3 )} [ l ~ f(tpi) ~ /0%J - f( £ P3 ) + f( e pi)f( e P2) 

+ f(e Pl )f(e P3 ) + f(s P2 )f(e P3 )l (88) 



where E = cu + e P2 +e P3 . After the substitution in eq. (|87| ) the desired form of the correlation 
function is found 

{f IH (e Pl )Ur = J^f(E)c(e Pl ,E), (89) 

where the third order quantum virial coefficient is defined as 

c(e pi , E) = 2 J2 I 1 - /( £ pJ - /( £ p 2 ) - /(e w ) + /(e P i)/(epa) + /(^i)/(^ 3 ) + /(£»)/(£«,) 



P2P3 



X 



(fc^ImT^lfc^^Re^E) + Imft(£) ^(A^lReT^)!*^) 



(90) 



and where 7£ = [12 — e pi — e P2 — e P3 + 177] is the three-particle resolvent. 
Collecting the results, the virial expansion for the density reads 

"tot = /|Sr W^ £ P ) /(^) + 6 (^> £ P ) 0(*O + c ( s > £ v) f( E )] ■ ( 91 ) 

This expression extends the quantum virial expansion for equilibrium Fermi-systems to terms 
of the third order in density. We note that the last term, as in the two-body case [|23|-|25| , can 
be separated further into bound and scattering contributions, which can then be expressed 
through the in-medium three-body scattering phase-shifts. 
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VI. SUMMARY 

In this paper, we derived coupled quantum kinetic equations for Fermi-systems in the 
framework of real-time Green's function theory. The quasiclassical double-time Kadanoff- 
Baym equation was reduced to two single-time Boltzmann-type equations for quasiparticles 
and off-mass-shell excitations which are coupled via generalized scattering integrals. The 
part of correlations related to the off-mass-shell excitations are accounted for via an expan- 



sion of the spectral function up to the next-to-leading order terms in spectral width [ JT,22 
which is regarded as a small parameter. Although the treatment of the drift part is concep- 
tually close to refs. pi| , |2"2"|] , we show that a different partition of the quasiclassical functions 



with respect to the orders of the spectral width can be exploited with the same efficiency as 
the separation in pole and off-pole contributions. In particular, the latter partition has the 
advantage that it fulfills the spectral sum rule at each order of the decomposition and, in 
addition, closely parallels the equilibrium treatments based on the same approximation to 
the spectral function p3|-p5|, thus providing a consistency check for the variables treated in 
the local equilibrium limit (e.g. scattering amplitudes). 

The collision integrals contain contributions from three-body scattering. The three- 
particle scattering amplitudes are determined by mapping the Faddeev decomposition on 
the time-contour. The resulting in-medium three-body equations contain the effects of the 
statistical suppression of the intermediate state three-particle propagation and quasiparticle 
renormalization; in the vacuum limit they reduce to the well-known Faddeev equations for 
the three-particle problem in free space ||26|| . The amplitudes derived in this manner differ 
from those considered by Bezzerides and DuBois |28| in that they resum the complete series 
in the particle-particle channel and are free of divergences associated with the non-Fredholm 
nature of the kernel of the three-body integral equation for the amplitudes. Our equations 
appear to be consistent with the time-local limit of the 2-particle-hole equations derived in 
the self-consistent random-phase approximation [p^7|p9| [3l[ . 



The decomposition of the scattering integrals in the on- and off-shell parts is accomplished 
in much the same manner as for the scattering amplitudes. The zeroth order term is the 
familiar Landau collision integral in the T- matrix approximation P|, P^JT7[] . Four further 
terms represent first order virial corrections to this collision integral, with one of the incoming 
or outgoing particles being in an off-mass-shell state. Consistent with keeping only the first 
order terms in the width, the transition probabilities enter the zeroth order term of the 
virial expansion with the full off-shell contribution, while the first order terms acquire only 
on-shell (zeroth order) contributions. A similar expansion holds for the kinetic equation 
for correlated elementary excitations. The scattering processes that contain an off-shell 
particle in the final state can further be expanded, putting all final states on the mass-shell: 
one then finds that these processes effectively correspond to successive two-body scattering 
events connected by an intermediate off-shell propagation. 

In the equilibrium limit the equation of state has the form of a virial expansion and is 
truncated at the level of the three-body correlation. An explicit expression for the third 
quantum virial coefficient, which is expressed via the three-body scattering amplitudes and 
the three particle resolvent, is obtained. The latter result extends previous discussions of the 



quantum Beth-Uhlenbeck formula for the density including the two-body correlation [p3|-|2~5 
to the three-body case. 
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Numerical studies of three-body correlations at finite temperatures and densities in the 
context of excited nuclear matter using the full three-body equations either in the AGS or 
Faddeev formulation are underway. The in-medium nucleon-deuteron cross-section has been 
calculated recently [Q , and the relaxation times for the nucleon-deuteron system have been 



estimated ||34|| . The Mott-dissociation of tritium in hot nuclear matter is now being studied. 
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